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A set of once subtracted dispersion relations with imposed crossing symmetry condition for the 
TTTT D- and _F-wave amplitudes is derived and analyzed. An example of numerical calculations in 
the effective two pion mass range from the threshold to 1.1 GeV is presented. It is shown that these 
new dispersion relations impose quite strong constraints on the analyzed tttt interactions and are 
very useful tools to test the tttt amplitudes. One of the goals of this work is to provide a complete 
set of equations required for easy use. Full analytical expressions are presented. Along with the well 
known dispersion relations successful in testing the tttt S- and P-wave amplitudes, those presented 
here for the D and F waves give a complete set of tools for analyzes of the tttt interactions. 
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I. INTRODUCTION 

The introduction of crossing symmetry conditions into 
dispersion relations for the tttt amplitudes was promoted 
by Roy in 1971 1]. Since then, now well-known the Roy's 
equations have been successfully used to test the tttt S- 
and P-wave amplitudes [3-[T2j| . These equations with two 
subtractions have been used for example to analyze the 
low energy tttt interaction parameters Q and to elimi- 
nate the long standing "up-down" ambiguity in scalar- 
isoscalar tttt wave amplitudes [3, @ . Quite recently in- 
terest in these equations has increased significantly due 
to series of works by Bern and Madrid groups j7HlO| . In 
these analyzes, authors have used, inter alia, the Roy's 
equations to construct the tttt S and P wave amplitudes 
fulfilling crossing symmetry for effective two pion mass 
m^TT from the threshold to over 1 GeV. To perform it they 
had to describe the tttt S-, P-, D-, F- and G-wave ampli- 
tudes usin g p henomenological parameterizations below 
1.42 GeV 01 or 2 GeV [S^ and Regge amplitudes up 
to several GeV. In such a way these analyzes delivered 
prescription for constructions of unitary tttt amplitudes 
for many partial waves in very wide energy ranges. To- 
gether with the outcome of new precise data near the tttt 
threshold [ll] these works led also to very accurate de- 
termination of the mass and width of the /o(600) (or a) 
resonance and of the threshold parameters [4, 11]. 

Very recently the tttt S and P wave amplitudes with 
crossing symmetry constraints have been analyzed using 
not only the Roy's equations but also once subtracted 
dispersion relations [ll|, [T^] . Due to one less subtraction 
these new equations, called GKPY, proved to be much 
more demanding than the Roy ones. Above around mTTTr 
= 450 MeV the GKPY equations have much smaller un- 
certainties so they impose stronger constraints on the tttt 
amplitudes. 

In the works mentioned above only the S and P wave 
amplitudes have been directly fitted to the Roy's or 
GKPY equations. The higher partial wave amplitudes 
have been fitted only indirectly via their relations with 
the 5* and P waves within the dispersive equations. 



In this paper are derived once subtracted dispersion 
relations with imposed crossing symmetry condition for 
D and F waves, hereafter called OSDRdf ■ Together 
with the Roy's and GKPY equations they will form a 
complementary set of dispersion relations which can be 
very useful in testing the tttt S-, P-, D- and i^-wave am- 
plitudes. 

This paper is organized as follows: in Section II gen- 
eral structure and analytical properties of the OSDRjyp 
are presented and discussed. Individual components of 
these equations are analyzed in detail. The third sec- 
tion contains an example of numerical calculations. Full 
derivation of OSDRjjf and analytical expressions for 
their components are given in Appendices A and B. Dis- 
cussion of results and summary are in Section IV. 

Through the text, partial waves with orbital momen- 
tum I and isospin / are denoted by WI or just by W if 
the isospin docs not need to be specified. The W can be 
S, P, D, F or G for / equal to 0, 1, 2, 3 or 4, respectively. 



II. DISPERSION RELATIONS 

Dispersion relations, for example the Roy's and GKPY 
ones, relate real parts of given partial wave amplitudes 
with sets of imaginary parts of other ones. General form 
of dispersion relations with one subtraction for the D- 
and P-wave amplitudes reads: 



J ds'^ii' ^')Ini //; (s') + 4 (s) 



(1) 



where s — rn^^. The first component — — 2.5aQ)/24 
is so called subtracting term ST^ . Two remaining parts 
are called the kernel KT^{s) and driving DT^ (s) terms. 

Below s' — s^(j2, partial wave amplitudes //(s) can 
be expressed by their relations with the experimental tttt 
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phase shifts S{ (s) and inelasticities rjf {s) 



where 



cr s 



2ia{s) 



Ami 



(2) 



(3) 



For brevity hereafter the real parts of amplitudes on the 
left side of Eq. ([IJ will be called "output amplitudes" 
and the imaginary parts of amplitudes on the right side 
- " input amplitudes" . 

The angular momentum V of the input amplitudes goes 
from to 3 (the 5*, P, D and F waves). As was shown in 
[T^ the input amplitude for the G wave {V = 4) is very 
small and has negligible influence on the output. Because 
of the Bose symmetry the sums I' + /' and I + / for input 
and output amplitudes respectively, must be even. 

As presented in Appendix A the subtracting terms 5T/ 
in once subtracted dispersion relations are constant and 
are determined by combinations of real parts of partial 
wave amplitudes at the tttt threshold. As can be seen 
from threshold expansion 



Refl{k) = k 
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where the pion momentum k — v (s/4 — m^), the only 
nonzero amplitudes at fc = are those for the 50 and 
S2 waves. Their values at the threshold are scatter- 
ing lengths flg and a§, respectively. In Appendix A it 
is shown that in the case of dispersion relations for the 
D and F partial waves the only nonzero combination of 
these two scattering lengths in the subtracting terms is 
for the F wave. It is due to nonzero integral (from to 
1) of the Legendre polynomial for this partial wave. The 
scattering lengths aj] and Oq can be treated as an input 
and may be fixed using for example ChPT predictions 
(see e.g. Q) or can be fitted to data and to theoretical 
constraints (see e.g. [lol - [l^ '). 

In the kernel terms KTg (s) the products of the input 
amplitudes and of the kernels K^J (s,s'), defined in Ap- 
pendix A, are integrated over s' from the tttt threshold to 
s^„^ i.e. up to the energy where the phenomenological 
parameterizations of the phase shifts and inelasticities in 
input amplitudes (see Eq. ^) are quite well known. In 
practice s^^,^ ^ 1.4 — 2 GeV (see e.g. analyzes of the 
Roy's [| and GKPY equations 12]). 

Above s' = s'^ax the input amplitudes are parameter- 
ized using the Regge formalism. Integrals in this s' range 
are grouped into the driving terms DTg (s). Contrary to 
the kernel terms, the products of the input amplitudes 
and corresponding kernels in DT^ (s) must be doubly in- 
tegrated - over s' and t. This is due to the t dependence 
of the Regge amplitudes. 

As shown in Appendix A, the only singularities in Eq. 
([T]) are those at s' = s in the diagonal kernel elements 



i.e for 1 = 1'. Therefore one has to take principal value 
there. 

As was already pointed out in [l^ while comparing the 
GKPY and Roy's equations for the S and P waves, an 
essential feature of the once subtracted dispersion rela- 
tions is their slower convergence than in twice subtracted 
dispersion relations. In the case of one subtraction, the 
integrands in the KT^ (s) and in DT^ (s) behave as 1/s'^ 
for s' — > 00 while in the Roy's equations as 1/s'^. Due 
to this difference the input amplitudes at higher energies 
in the former equations enter with higher weights than 
in the latter ones. It is especially important for higher 
partial waves which contribute mainly above 1 GeV and 
therefore are undervalued in relations with two subtrac- 
tions. Since the output amplitudes for given / partial 
wave depend also on the input ones with V = / this ar- 
gument is reinforced in this analysis dealing with D and 
F output partial wave amplitudes. 

Constant value of the ST/ in OSDRdf leads to con- 
stant value of its errors which, one can expect, are much 
smaller (for higher s) than these which would be in anal- 
ogous relations with two subtractions. It is due to the 
fact that subtracting terms in these relations are not con- 
stant but are linear functions of s which leads to increase 
of their uncertainties with energy. Detailed comparison 
of the errors in the once and twice subtracted dispersion 
relations for the S and P waves (the GKPY and Roy's 
equations, respectively) can be found in 12j . 

The application range of the OSDRdf , Eq- (H]), is the 
same as of the GKPY equations in and comes about 
1.1 GeV. One can expect that below this energy only the 
DO output amplitude will show a clear increase with the 
energy due to presence of the /2(1270) resonance. The 
absence of any resonance in the D2 wave and the rela- 
tively large mass of meson p3(1690) in the Fl one, lead 
to rather small variation of amplitudes in these waves 
below 1.1 GeV. 



Wave 


TE 


ST/ 




OT/(s) 


DO 








ao + cqPo 


70 + (IqSq 


Fl 









5 + 71 + 5i 


D2 








02 + C2/32 


72 + d2<52 



TABLE I: Comparison of the threshold expansion TE (Eq. 
(HI) results for s — >■ 4m^ and of the threshold behavior of 
subtracting S'T/, kernel KTI{s) and driving DTl{s) terms 
for D and F partial waves. The ai and 7/ are terms of the 
order 0{s — Am%) and 5i of the order 0(s — 4mJ)^. Values 
of the A, B, ci and di constants are explained in the text. 



In practical applications of the OSDRjjp it is very in- 
teresting and useful to compare the threshold behavior 
of the S'T/, KT^^is) and DTl{s) for D and F waves. Ac- 
cording to the threshold expansion the sum of these 
components should vanish at s = 4mJ for all these waves. 
Comparison of the threshold expansions of the subtract- 
ing, kernel and driving terms is presented in Table HI For 
the DO and D2 waves the zero order parts in all compo- 
nents are equal to zero, which automatically ensures cor- 
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rect behavior of the full output amplitudes at s = 4m^ . 
Of course sum of the first order parts aj and 7/ must 
be equal to zero and sum of the c/ and dj should give 
corresponding scattering lengths. In case of the .Fl wave 
the nonzero value of subtracting term must be canceled 
by sum of nonzero values of KT^ {Arn^) and DT^ {Am^). 
Sums of the first and second order parts should give zero. 
It is worthy to emphasize here that such cancellations 



demand strong and proper mutual relations between the 
amplitudes for all waves integrated over very wide energy 
range m the KTl{s) and DTl{s). These relations for 
the zero order parts of the F\ wave, i.e. for A and B 
in Table |T1 involve also the threshold parameters of the 
lowest partial waves (5*0 and 5*2) and are expressed by 
the well known Olsson sum rule [U] 



- lal) 



ds' 



4771?. 



2/m/^(s') + 10/m/2"(s') +9/m//(s') + 2llmf^{s') - Mmf^{s') - 2Mmfi{s') 
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s'(s' - Ami) 



(5) 



where Fl{s',t) is the isospin 1 amplitude in the t channel (see Appendix A). 



Of course an identical relation may also be obtained from 
the GKPY equations. In practice this sum rule gives a 
chance to verify an accuracy of used parameterizations 
of many input partial wave amplitudes in a very wide 
energy range and validates the choice (or result of a fit) 
of scattering lengths ag and Oq. 



III. NUMERICAL ANALYSIS 

As an example of practical application of the 
OSDRdf, the D and F output amplitudes have been 
calculated using the input ones from [12J. These input 
amplitudes were obtained there in dispersive data anal- 
ysis using forward dispersion relations, sum rules, the 
Roy's and GKPY equations. The D and F waves have 
not been directly fitted to any dispersion relation. 

Figures [1] to [3] present niT„T distributions of the output 
amplitudes and the real parts of the amplitudes whose 
imaginary parts have been used as the inputs in Eq. ([T|) . 
The error bands of the output amplitudes have been cal- 
culated using Monte Carlo method for all 53 parame- 
ters used in |12| to parameterize the low and high energy 
behavior of all input amplitudes. Assuming Gaussian 
distributions of these parameters they have been varied 
randomly 10'^ times within their left and right 3cr which 
resulted in three Gaussian distributions for the three out- 
put amplitudes. This procedure has been carried out for 
25 values of s between the tttt threshold and 1110 MeV. 
The errors of the output amplitudes have been deter- 
mined by widths of the Gaussian functions fitted to the 
left and right sides of these distributions at each s in- 
dependently. Of course, due to the large number of pa- 
rameters, these errors are very similar to those obtained 
as square root of sum of squares of individual errors of 



each parameter. However, the Monte Carlo method has 
been chosen to present possible asymmetries of final er- 
rors caused by correlations between varied parameters. 

As is seen on the figures [T] to [3] the error bands of 
the output DO and D2 wave amplitudes go to zero for 
s Arn^. The nonzero errors of the Fl wave at the 
threshold are due to the nonzero value of subtracting 
term and in fact represent its error. In Table U this 
constant ST^ is a linear combination of two scattering 
lengths of the 5*0 and 52 waves. Therefore, the full er- 
ror of the output Fl wave amplitude at the threshold is 
completely determined by the uncertainties of these two 
threshold parameters. Convergence of the DO and D2 
output amplitudes to for s — 4mJ is well seen. In 
the case of the Fl wave, the output amplitude goes to 
(— 1 ± 7) • 10"'' which is compatible with the predicted 
zero value in Table HI The error has been calculated in 
the same way as the errors of the output amplitudes. 

As the D and F wave amplitudes have not been di- 
rectly fitted to any dispersion relation in [13], the fact 
that their input and output amplitudes for m^^ > 
800 MeV differ by more than one a is not surprising. 
One can expect, however, that the use of the presented 
here dispersion relations in more complete analysis of 
the TTTT amplitudes e.g. in the analysis of the Bern or 
Madrid group would decrease this difference significantly. 
It would also increase the weight of the theoretical con- 
straints imposed on these and on other wave amplitudes, 
and therefore would diminish their uncertainties. 

Vastly different scales in the figures for the DO and 
for D2 and Fl waves reflect significant differences in the 
sizes of the amplitudes. As was explained in Section |TT] 
these differences are due to the presence of the meson 
72(1270) in the DO wave close to the studied s region 
and lack of such states in the D2 and Fl waves. 
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FIG. 1: Input (solid line) and output (dashed line) for the 
DO wave amplitude in Eq. ((TJ. Gray band represents errors 
of the output. 
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FIG. 2: As in Fig. [T]but for Fl wave amplitude. 
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FIG. 3: As in Fig. [T]but for D2 wave amplitude. 
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FIG. 4: Components of the output DO wave amplitude: kernel 
term (solid line) and driving term (dashed line). Gray bands 
represent their errors. 
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FIG. 5: As in Fig. |4] but for the F\ amplitude. Horizontal 
dashed-dotted line represents subtracting term. 
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FIG. 6: As in Fig. |4]but for the _D2 wave amplitude. 



Figures. |4] to |6] present the ttItttt dependences of the 
ST/, Xr/(s) and DTl^s) components for the D and 
F wave amplitudes. Error bands have been calculated 
for each term separately as for the output amplitudes in 
Figs. [l]to[3l Again, due to the presence of the /2(1270) 
resonance much larger values of the i^r/(s) and DT/(s) 
for the DO wave than for the _D2 and FX ones are visible 
(note different scales on the figures). 

For all three waves the errors of the kernel terms are 
significantly smaller than those of the driving ones. This 
is due to the fact that the input amplitudes in KTl{s) are 
more precisely known than those in DT^{s). The errors 
of the latters are determined by errors of the input am- 
plitudes above s'^^^. These amplitudes are, however, less 
precisely known than those at lower energies. In particu- 
lar, small errors of the kernel for the DO wave amplitude 
are due to presence of the well known resonance /2(1270). 

In Fig. [S]for the Fl wave, the straight line represents 
nonzero subtracting term value of — 0.01371o!qoo9- ^^s 
errors are completely determined by the uncertainties of 
the scattering lengths of the 50 and 5*2 waves. If the 
amplitudes of these waves at the threshold are not fixed 
but fitted (as it is e.g. in [l^]) then of course, these un- 
certainties become functions of the errors of amplitude 
parameters. According to what was shown in Table U 
the kernel and driving terms also have nonzero values 
at the threshold to compensate the subtracting term. 
Threshold values of the kernel and driving terms are: 
A = 0.0113 ± 0.0002 and B = 0.0023 ± 0.0002, respec- 
tively. 



IV. CONCLUSIONS 

A set of once subtracted dispersion relations with im- 
posed crossing symmetry condition for the tttt D- and 
_F-wave amplitudes has been derived and analyzed. An- 
alytical structure of these equations and of their compo- 
nents have been studied and described in detail. It was 
shown that integrals in these equations converge slower 
than integrals in the twice subtracted dispersion relations 
(e.g. in the Roy's ones for the S and P partial waves). 
Thanks to this, the input amplitudes are less suppressed 
at higher energies which is essential for higher partial 
waves becoming significant only above 1 GeV. One can 
expect that use of the one subtraction will lead to the 
generation of smaller uncertainties of the output ampli- 
tudes than those created using two subtractions. In the 
latter case, the errors of the subtracting terms, being not 
constant but linear functions of s, grow with s. There- 
fore, the one subtracted dispersion relations provide more 
demanding tests for the tttt amplitudes. 

Apart of the derivation of these equations, a first prac- 
tical application has been presented. The input ampli- 
tudes in the low and high energy region have been taken 
from fit to experimental data and from direct, for the 
S and P wave amplitudes, and indirect fit to dispersive 
theoretical constraints in [l2| . It has been shown that be- 
cause of the D- and i^-wave amplitudes were not fitted 
directly to any theoretical constraints there, their inputs 
and outputs, calculated in this paper, differ sizable above 
about 800 MeV. One can expect, however, that the com- 
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patibility between them will be significantly improved if 
presented here dispersion relations are used in fits of the 
TTTT amplitudes. 

It was shown that nonzero value of the subtracting 
term in the Fl output amplitude creates opportunity to 
relate low and high energy behavior of many partial wave 
amplitudes with the threshold parameters of the lowest 
ones. Therefore one can expect that these equations may 
help to decrease uncertainties of the tttt amplitudes and 
errors of these threshold parameters. One can also ex- 
pect that the theoretical constraints imposed by these 
equations directly on the D and F wave amplitudes and 
indirectly on lower ones can lead to slightly more precise 
determination of the /o(600) and /o(980) parameters. 

In conclusion the above derived and analyzed set of 
equations is easy to use and can be very helpful in future 
analyzes of the tttt interactions. Together with disper- 
sion relations for the S and P waves - the Roy's and 
GKPY equations with two and one subtraction, respec- 
tively, they can be used as complementary set of equa- 
tions to test TTTT amplitudes from the threshold to about 
1.1 GeV. 
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Appendix A: Derivation of the once subtracted 
dispersion relations 

Let us define following isovector for the scattering am- 
plitudes F^(s, t) of isospin / 



F(5,t) 



F^{s,t) 



(Al) 



Then a once subtracted dispersion relations can be ex- 
pressed by 



Re F{s,t) = Re F(so,t) 



So 



(A2) 



f ^ , lmF{s',t) ^ f ^, 



Im F{s',t) 



(s'-so)(s'-s) 



where sq is the subtraction point to be defined later. 

The first and second integral in Eq. ()A2p are taken 
on the real s axis along the right and left hand cuts of 
F^{s,t), respectively. For any s values along these cuts, 
one should take the principal values for these integrals 
(see e.g. the discussion in Appendix B of 12]). 

Performing the substitution 



and using the crossing symmetry relation 
F{u',t)^CsuF{.s',t), 



(A4) 



the left hand cut integrals can be recast in terms of the 
right hand cut ones 



Re F{s,t) = Re F{so,t) 



So 



(A5) 



X + ds' 

4m?. 



lmF{s',t) lmF{s',t) 



(s'-so)(s'-s) 



' (s' - uo){s' - u) 



with u — 4m^ — s — t, uq — 4m^ — sq — t and with the 
crossing symmetry matrix Csu defined in Eq. (|A8|) . 

Using the s < — > t crossing symmetry relation one can 
express the subtracting terms F^{so,t) by 



F(so,0 = CstF{t,so). 



(A6) 



Following now the same procedure which was used in 
derivation of the KeF^{s,t) in Eq. (|A5|) . one can get 
similar, once subtracted dispersion relations for F^ {t, sq)- 



t-to 



Re Fit, So) ^ Re F(io,so) 

TT 

oo 

X 

Am' 



(A7) 



OO 



Im F{s',so) 
is' ~t„){s' -t) 



Im F{s',So) 
(s' - uoo){s' - uo) 




4m: 



(A3) 



with Moo ~ 'itTi'l- — Sq — to and to = 4mJ — so — u. The 
crossing matrices Csu and Cst used in Eqs. (jA4l) to (|A7[) 
read 



(A8) 



It is worthy to mention here that the set of equations 
(IA7|) can be easily obtained from (lASp using the following 
replacements: s t, t so, so ^ to, u uo and 
Mo Woo- 

The problem of the convergence of integrals (|A5I) and 
(|A7P has been analyzed and described in [l2|. It was 
shown there that, due to the Pomeron contribution from 
the /( = channel, the integrals along the left and right 
hand cuts are divergent when taken separately. How- 
ever, taking both cuts into account simultaneously, a can- 
cellation occurs and the integrands decay as 1/s'^ when 
s' — > oo which ensures convergence of the integrals. In 
the Roy's equations with two subtractions corresponding 
integrands behave like 1/s'^ which ensures convergence 
of each integral separately and of course faster, than in 
case of once subtracted equations, convergence of their 
sum. 

Substituting now (jA6p and (jA7l) into (jA5l) and follow- 
ing the Roy's original choice: to = 4toJ and so = one 
gets 
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Re F{s,t) = ujCsta 

t - Ami 



ds'Cs 



4mJ 



ImF(s',0) Cs Jm F(s',0) 



{s' -t){s' - Ami) s'(s'-Mo) 



+ 



Im F{s',t) Csul-oi F{s',t) 



s'{s' — s) (s' — uq){s' — u) 



where to is some constant (for example 32tt in ^ and S/tt 
in [i3) and 



^de/ F(4m^,0) 
a = 



"0 



(AlO) 



is a vector with, the elements of which are scattering 
lengths of the SO and iS'2 wave amplitudes and are defined 
by the threshold expansions given in Eq. Q. 

Projection of the vector F{s, t) on partial waves //(s) 
is given by 



dxPi{x)F{s,t), 



(All) 



r 



(|A14I) and can be shortly written as 
R-e fi,{s) = Csta 



4ml 



(A9) 



(A15) 



where = J dxPi{x) and Re ff> (s) is the hight en- 



ergy part defined later. For simplicity all integrals in 
this equation are grouped in 

Ku'is,s) = (2f + 1) (A16) 
X I Ku,{s,s')i - Lii,{s,s')Csu 



where 



Ms)^ fl(s) 
ffis) 

Pi{x) are Legendre polynomials and 



(A12) 



t = 



(s - Aml){x - 1) 



(A13) 



Here, due to the symmetry of the integrands, the integra- 
tion limit was taken from to 1 instead of from -1 to 1. 

In the projection of the isospin amplitudes on partial 
waves it is convenient to split integrals in Eq. (IA9p into 
low (for s < s'^^^) end high energy part (s > s'^^J. 
Value of s^^j, is determined by limited range of appli- 
cability of phenomenological parameterizations for the 
input amplitudes. Typically is s'^^^^ ^ 1.4 — 2 GeV (see 
for example and |12|). 

In the low energy part the input amplitudes on the 
right hand side of (jA9[) are expressed by 

Fis,t)^LuJ2(^l + ^)Pii^)fii^) (A14) 

Analytical expressions for the output partial wave am- 
plitudes are obtained from the relations (IA9|) . (jAllI) and 



The Ka,{s,s'), Lwis.s'), Mi{s,s') and Ni{s,s') 
kernels read 



TTS'iS — S' 



Lu'{s,s')) = -UxPiix) 

TT J U (U — S) 





j dxPi{x)Pi, {y) , (A17) 
(A18) 



Mi{s,s') 
and 



1 } f — 4r77^ 

, , , / dxP,ixf-^{Al9) 
n(s — 4r7?4) J s' — t 



i 

Ni{s,s') = idxPiix)^^ 
its' J w 



(A20) 



where y — [u' ~ t)/{u' + t). Integrands of the kernels 
Lii'{s, s') and -/V;(s, s') have singularities at u' = which 
for s' < ■^{s + Am'^) is in the range of integration. There- 
fore principal values are taken for these integrals. 

The full analytical expressions for the K^^ elements 
of Kiii{s, s') are presented in Appendix B. 
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In the high energy parts of the integrals in Eqs. (jA9l) . 
the input amphtudes can be expressed by t-channel ones 
Ft{s,t) using the Regge parameterizations. The relation 
between amplitudes in the s and t channels is 

Fis,t)^CstFt{s,t). (A21) 



Following it, one can conclude that high energy partial 
wave amplitudes derived directly from Eqs. (|A9p read 



oo ^ i 

Re fi'^'is) = j ds'l^ J dxPi{x) 



C.S 



s'(s' — s) (s' — uo){s' -~ u) 
(^Mi{s',s)i-Ni{s',s)Csu) CstImFtis\0) 



CstImFt{s',t) (A22) 



Finally, the full expression for given output partial 
wave amplitude can be written as in Eq. ([T]): 

i?e//(s) ST/ + KTf'is) + DT^^s) (A23) 

where for the subtracting term: 

5T/ = Cstct = -^(a[S - Ial)dn6i3, (A24) 
and for the kernel terms: 

The driving terms DT'{s) are given by Eq. (|A22|) . 

Appendix B: Analytical expression for kernels 

Following work [15] for the Roy's equations, one can 
express eighteen kernels K^^, (s, s') in Eq. (jA25P as func- 
tions of the four kernels Kii'(s,s') and eight combina- 
tions of the L;//(s,s'), Mi{s,s') and Ni{s,s') ones from 
Eqs. ((ITtT) to ^M0^ 



where 



1^00 

■"■21 
1^02 
■"22 


~ K20 — ho/i, 
= 3/21, 

25 J 


^2°0^ 
K'2^ 
■"23 


- -3^20 
= HK22 
= 7/23, 


- K22), 


■"30 
■"31 
■"32 


= /30/3, 

- 3(/f3i - 5/31), 

25 T 


■"30 

■"32 
■"33 


- -6-'30 
5 r 

- 3^32, 

= 7(^33 


- 1-^33), 


■"20 

7>'21 

■"21 
1^22 
■"22 


= -/20/3, 
3 T 

- -2-'21, 

= 5(if22 - ^/22), 


■"20 
1^20 
■"22 
1^21 
■"23 


= ■?^20- 

— — 3^22 

- ^2^3 


/20/6, 



^20 


= ^^20 


-M2 


+ N2, 


/2I 


= ^^21 


+ M2 


~N2, 


^22 


= -^22 


-M2 


+ N2, 


^23 


= -^23 


+ M2 


~N2. 


^30 


— L30 


+ M3 


+ N3, 


/3I 


= i31 


-Ms 




■^32 


= ^32 


+ M3 




■^33 


— 


-Ah 


~N3. 



(Bl) 



Full analytical forms of the kernels Kni are 



i^20 

K22 
K31 
K33 



0, 

(s - Aml)s{7s - 15s' + 32m^) 

407r(s -s')(s' -4to2)2s' ' 
s 



87rs'(477l2 - s')' 

s 



(B2) 



567rs'(4m2 — s')^{s — s') 
[-512m^ + 8s^ - 64m^(s - 7s') - 35s2s' 
42ss'2 - 7s'3 -I- (44s2 - 56ss' - 84s 



/2^ 



The other (even nonzero) Ku' elements are not given 
here as they are multiplied, in (jA25[) . by zero //(s) am- 
plitudes. It is caused by the Bose symmetry for two pions 
according to which the sum / + I must be even. It also 
leads to the same condition for sum I + I' in the equa- 
tion (|A15|) . The diagonal kernels K22{s, s') and i^^33(s, s') 
contain singularity at s = s' which is the only type of sin- 
gularities in presented equations. 

In order to simplify the analytical forms of the /;/' el- 
ements, some terms which appear at least twice or help 
to reduce these formulas are given below: 
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ai = s-Aml, a[ = s'-Aml, a2 = 2s' - Sm^, 03 = s + 2s' - 4m^, 04 = 2s' - s - 4m^, 
s' \ ^ , / 2s' \ ^ , / s + 2s' - 4m2 \ 



Li = log 

L4 = log 



S + S' — 4777,2 

2.s' - 8;;4 
-s + 2s' - 4ml 



L2 = log 



s + 2s' - 4ml 

L5 = log (2(-4m2 + s + s')) , 



, L3 = log 



'^V2(s + s'-4m2)j 



(B3) 



/i = 16m^ + s2 + 8m2(2s-3s')-6ss' + 6s'^ 

/2 = 16mt + s^ + 6ss' + 6s'2-8m^(s + 3s'), 

/s = 64777^ + + 48m^(3s - 4s') - 12s2s' + 30ss'2 - 20s'^ + 12to^(3s2 - 12ss' + lOs'^), 

/4 = 64777^ - - 12s^s' - 30ss'2 - 20s'3 - 48777,t(s + 4s') + 12777^ (s^ + 8ss' + lOs'2), 

/s = 192777,^ - 3s3 - 148s2s' + BOOss'^ - 480s'3 - 16777^(9s + 28s') + 4m^(9s2 - 304ss' + 360s'2), 

fe = -64ml + 48(s + 4s')m^ + s^ + 20s'^ + 30ss'^ - 12 (s^ + 8s's + lOs'^) ml + 12sls', 



(B4) 



f - "1 



9280m^ - 16(287s + 748 s')m^ + 4 (l39s2 + 896s's + 1080s'^) to^ + 3 s^ - 480s'^ - 600ss'^ - 148s2s' 



+ 48/4^2 + 

s' 

The analytical expressions of the Iw elements can then be expressed as 

ho = - o ^3 (- 3aia3 - /2ii), 
3ai''7r ^ ' 



(B5) 



'21 



ai^a'iTT 



64ml - 16(4s + 7s')rnt + 4 (ss^ + 20s' s + 12s'^) ml - 2s^ - 6s'^ - ISss'^ - IBs"^ s' 



+ 9sai^ + 24a[mlai - 12s^ai + 72m^s ai - 12a[s'ai - ISss'ai - 2ai/2i2|, 



(B6) 



122 



5 f 3oi {12ml - 3s - 4s') + 3 ai {4ml + 3s - 4s') + 4/1X4 - 4/2^2 



6tt 1 



+ 



2s 

~T2 

a'l 



ai^s 



9g/ 6 (32m^ - 8(s + 4s')m2 + + 7 s' V 6ss') 

^a'l^ /1L4 + is ^256m* - 512 (s + s')m% + 16 (Ws^ + 52s's + 19s'^) m% 



(B7) 



- 8 (9s^ + 43s's2 + 43s'^s + 9s"'') + 6 s^ + 6s"' + 42ss'^ + 73s2s'^ + 42s^ ) ^ 



-^23 = 



4a\^'K 



6ai {4ml + 3s - 4s') 



5824m^ - 48(79s + 168s')mt + 12 (95s^ + 348s' s + 240s'^) ml - 115s^ - 312s'^ - 792ss 
1024mi° - 256(14s + 9 s')m\ + 64 {h^s^ + 108s's + 27s'^) m% 

"1 

16 (92s^ + 333s's^ + 252s'^s + 37s'^) + 4 (70s^ + 384s's^ + 531s'^s^ + 236s'^s + 24 s'*) ml 



'' 660 s^s' 



+ 773^3 



- 20s^ - 6s'^ - 78ss''' - 253s2 s'^ - 312s3s'^ - ISOs^'s' 

- 2 (ai (-36 + 9s + 12s') + 4/2X2) >, 



(B8) 
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72 ai^a'iTrs 



— ^< 120%ss'ai^ 



a[ ^192 ml - 16(9s + 285')^^ + 4 (ds^ + 176s's + 360 s'^) - 3s^ 
- 480s'^ - 360ss'^ - 1485^ s'^ + s' |^9280to^ - 16(287s + 7485')™^ + 4 (l39s2 + 896s's + 1080s'^) ml 



+ 3s^ - 480 s'^ - 600ss'^ - 1483^5' 



/2 



ai - feA8a[s'Li }, 



(B9) 



731 = — — i- — -\ log(s')f - 3072a; s'm^ + 9216a'i s'^m^ + 2304a'iss'm^ - 5760a[s'^ ml - 4608a'iss'^m^ 
- blQa'^s^s'ml + 960 o'ls'" + 1440a'iss'^ + 576ais^ s'^ + 48a'iS^s') + 48a'ilog(2)s76 



+ ai 



s' ^ - 9280m^ + 16(863s + 748s')m^ + 4 (^-331 + 240ais - 3296s's - 1080s'^) ml + 93s^ 



+ 480 s'^ + 2520ss'^ + 80ai^s + 11085^5' - 120 axs{s + 3s') 



(BIO) 



+ a\y- 192m^ + 16(9s + 28 s')m1 - 4 {^s^ - 304s's + 360s'^) + 3 + 480s'^ - 600ss'^ + 148s^s' 

256m^ - 64(5s + 13s')m^ + 48 (35^ + 17 s's + 14s'^) - 4 (75^ + 63s's2 + 114 s'^s + 50s'^) ml 



- 48 s' 



+ 2s^ + 20s'* + 70s s'^ + 72s^s'^ + 25s^s' 



is - 48 (4m^ + a'l - 2s - s') s'/g log(a3) \, 



I32 



5 . _ 4^ 
72ai47rV^' a'" 



1024mi° - 256(9s + 14 s')m^ + 64 (^27s^ + 108s's + 55s'^) 



(Bll) 



- 16 (37s3 + 252s's2 + 333s'^s + 92s'^) + 4 (24s'' + 236s's3 + 531s'^s2 + 384s'^s + 70 s'"^) ml 

- 6s^ - 20s'^ - ISOss'" - 312s^ s'^ - 253s^s'^ - 78s*s' 



6ais 



5824 ml - 16(73s + 668s')m^ + 4 (93s2 + 736s's + 1220 s'^) ml - 31s^ - 640s'^ - 820ss'^ - 388s2 s' | 



and 



733 = ^ /,-^ 



96768 ml - 192(277s + 1177s')m^ + 16 (l513s^ + 8517s's + 9348 s'^) (B12) 
4 (^1145s^ + 10221s's^ + 17076 s'^s + 9540s'^) ml + 325s^ + 3360s'* + 9300s s'^ + 10368s^s'^ + 4035s^ s' 

4096m^2 - 15360 (s + s')mi° + 1445s^s'^ + 972s^s'^ + 270s^s' + 768 (23s^ + 67s's + 23s'^) ml 

- 192 (49s^ + 257s's2 + 257s'^s + 49 s'^) ml + 144 (iSs'' + Uls's^ + 253s'^ + Uls'^s + 18s'*) 

- 36 (lOs^ + 106 s's^ + 283s'^s^ + 283s'%^ + 106s'*s + 10 s'^) ml + 20s^ + 20s''' + 270ss'^ + 972s2 s'* | . 



48 
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